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© IMAGE PROCESSING APPARATUS. 

© The image processing apparatus prepares calculated image of high quality by effecting calculation in 
reduced amounts at the time of y-fitting within short periods of time maintaining high precision. The y-function j S 
transformed into a linear equation through logarithmic transformation in order to decrease the amount of 
calculation for y-fitting. the measured data after the logarithmic transformation are weighted, and disturbance is 
removed from the measured data to improve precision for Y-fitting. The calculated image data are operated 
based on a y-function curve fitted for each of the pixels, and the data are interpolated for those pixels that could 
not be calculated. jn order to obtain smooth calculated image. 
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The present invention relates to an image processor in which gamma function curves are set which 
approx.mate to curves showing the movement of a contrast medium on a tomogram image taken by a 
tomograph and various characteristics of the set gamma function curves are constituted as a computer 



(Background Art) 

A tomograph is a device for obtaining a tomogram image of a subject body based on X-rav 
transm.ss.on data obtained by irradiating the subject from its surroundings with X rays. As one oVtte 
techn,ques of. tomography there is a method for diagnosis of a morbid part as shown n the follow^ a 
water-soluble contrast medium is injected into a vein and an image of a blood vessel or an internal^ 
e whch changes with the lapse of time are observed. The time serial measurement data of a cTrtl po^ni 
or a part of concern in a tomogram image are said to form a gamma function curve as shown in Fiq 3 

o W l n a L 7h *? ^ 88 3X18 °' 3bSCiSSaS and ,he Change if1 a CT numbar * taken aT^ea^s of 
ordmates. therefore .t is necessary to set a gamma function which approximates most precisely to the 

S M IS T T l ° n t0m ° 9ram imaQeS 0btained in respective *n" <°" ^spective 

slices), and tiiey are connected w,th a gamma function curve 21 which approximates most precisely to the 
time senal change of the measurement data. The gamma function is expressed in equation (1). 

f=K(t-t 0 )- exp \-B (t-to ) } -(I) 

As a method for gamma fitting, there is a method in which parameters for making an error minimum is 

bought ,nto a norma equation by transforming the gamma function into another function, for example, a 
logarithmic function. In the former case, that is. in an asymptotic method a great many calculations are 

m e *«H T TT r - th ? S ° MOn " 0t be 3,WayS ° bteined - ■» ,atter «»■ ■« a transformation 
T^J?^ 9 r T ' H n ?° n ' th6re t0 ° ' ar9e err0fS f0r pracBcal use in the area ^ va.ues of data 

S^Tnni" ™ X? T 1 S * UU t0 Perf ° rm aCCUrate 9amma fiWn9 with a smaN amount of calculations, and 
•t was not possible to obtain a high quality computer image using the gamma function curves. 

(Disclosure of the Invention) 

n am t n a °« b i eC! ? 6 PreS6nt inVenti0 " iS t0 rea ' i2e an image Pressor in which an amount of calculation for 

? m T ' S r reaSed ' 3 9amma Mn9 ° f hi9h precision is performed in a *<« time and a computer 

image of high quality can be formed. 

..f -^-.T, 396 P , f0Ce$$0r a l CO fd ' n9 t0 "» PreSent inventi0n functions as sh own in thfi f """ wi "0- ^ amount 
IS^-H^f^ 9 iS deCTeaSed by lrans,ormin 9 a gamma function to a first degree equation 

SlT^lS^^f 0 ^ 5 ? 3CCUraCy ° f 9amma f ' ttin9 i$ impr ° Ved by Anting measurement 
smLh Tol L ,0 9 ar '«™ c unction, and by removing disturbance in the measurement data; a 

whS J! / ma " Ca " ^ f0fmed by ° pera,ing COmpu,er image da,a ,ram gamma function curves 
riemente ^ P ' CtUre e ' ementS ' 30(1 b * interp0,atin 9 the ^ata for incalculable picture 



(Brief Description of the Drawings) 

Rg. I is a conceptual block diagram of an embodiment of a gamma fitting device to be used for an 
•mage processor according to the present invention; 9 ^ ° 3 

Rg. 2 is a conceptual block diagram of an embodiment of an image processor according to the present 
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invention; 

Fig. 3 is a graph showing a gamma function curve which is approximate to measurement data; 

Fig. 4 is a graph showing an example of input data to a gamma fitting device; 

Fig. 5 is an illustrative representation of coordinate transformation by a data clip device; 

Fig. 6 is a graph showing a relational curve between a real space and a logarithmic space; 

Fig. 7 is a graph showing input data containing disturbance; and 

Fig. 8 is an illustrative representation of a method for deciding to. 



io (Most Preferred Configuration for Executing the Invention) 

In Fig. 1, 1 is a data clip device which replaces a datum having a negative CT number to that having a 
positive number close to zero when the time serial measurement data of a certain point or a part of concern 
on a plurality of- tomogram images obtained by slicing in a plurality of times of scans are transformed to the 

is data in a coordinate system which is based on the time and the CT number of the leading data. The data 
clip device processes input data as shown below. An example of input data is shown in Fig. 4, and the 
process for the data will be explained. Each point in the figure is plotted based on the data obtained by 
slicing at each point of time. The numbers attached to the data are data numbers, and data 1 express the 
data of CT number CT1 obtained at a time Tl in the first slicing, and data i express the data of CT number 

so CTi obtained at a time Ti in the No. i slicing. The data clip 1 transforms the data shown in Fig. 4 for gamma 
fitting to the data in a coordinate system based on the time Tl and the CT number CTI of data 1, a leading 
data, that is, the data just before the arrival of a contrast medium. The coordinate-transformed data are 
shown in Fig. 5. In the graph shown in Fig. 5, the increment of the CT number of data 2 is negative, so that 
the data clip device 1 replaces the negative value with a certain value A close to zero which is set 

25 beforehand (A > 0). The coordinate transformation is. therefore performed as shown in equation (2). 

(T1 , CT1 ) ~> (t1 , M1') — (2) 

where t1 = T1 -Tl 

Ml' = CT1 - CTI (when CTI - CT1 1 A) 

= A (when CT1 - CT1 < A) 
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A part 2 is a logarithmic transformation device which sustains an expression of a gamma function in 
equation (1) in the form of a first degree equation being transformed to a logarithmic function, and 
transforms the output data of the data clip device 1 to a logarithmic function and substitutes them 
successively for the coefficients in equation (1). The device sustains the expression of gamma function in 
equation (1) in the form of a three element first degree equation being transformed to a logarithmic function 
as shown in equation (3). 

M=ax+byfcz (3) 
where M— log « f 
x=!og * K 
y = a 
z = & 
a=l 

b=log « (t — to ) (Hereinafter abbreviate lqge as log) 
C=~ (t-to ) 
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In the data clip device 1. the respective values M1. bl and cl in equation (3) are given by the 
coordinate-transformed output data (tl. MV). Therefore, equation (4) is obtained from equation (3). A proper 
value is applied to to until the value of it is decided later. 

x-f-bi y-f-ci z— Mi =0 '"(4) 



A part 3 is a weighting device to make the accuracy of gamma fitting in a logarithmic space equivalent 
to the accuracy of gamma fitting in a real space. Even if gamma fitting is performed by the output data of 
the loganthmic transformation device 2. it is a phenomenon in the logarithmic space, and the fittinq is 
largely influenced by the measurement data of small values, therefore the accuracy for approximation in the 
real space is degraded. The object of the weighting device 3 is to improve the above-mentioned problem 
Let ei be the residual at a measurement point i, when the left side of equation (4) does not become zero In 
general for gamma fitting the values of x. y and z which make the square value sum of ei's minimum are to 
'5 be found as shown in equation (5). 

n n 

2 ei ' 4 = 2 (x-hb, y-t-ci z-M, ) 2 ... (5) 
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Since the evaluation of equation (5) is made concerning a logarithmic space, in the present invention to 
jmprove the fitting accuracy in a real space, equation (5) is evaluated by a weighting function wl. a function 

2 e t 2 = 2 (w, e, ' ) 2 
i«l H 

n 

= 2 { (x-t-bi yrci z-M t ) w, } 2 — (6) 

35 t ' 0t ° f ,^ in9S Ca " 56 considered about the above-mentioned weighting function wl. but an example of 
them w,ll be explained referring to Fig. 6. In comparison with a data set X in a real space, consideration is 
made on the log X which is obtained by logarithmic transformation of the data set X Fig 6 is a qraoh 
showing a curve of Y = log X corresponding to the data X in a real space. In the figure, let (X0. logXO) be 
the coordinates of P. then the change sensitivity in the ordinate Y for the change in the abscissa X of data P 

40 can be expressed by a differential coefficient 1/X0 at point P on the curve log X. Therefore it is made 
possible to obtam gamma fitting which is not influenced by nonlinearity generated in logarithmic" transforat- 
ion for an arbitrary data X by multiplying the reciprocal of a differential coefficient by the residual in a 
oganthmic space. Because of this, in the case of transformation to a logarithmic space, the weighting 
function w (X) can be decided to be w (X) = l/(l/X) = X. Let w, in equation (6) be wl = MV and let the 

J5 expression be sustained in the weighting device 3. and by substituting input data t1. Ml' successively the 

_ . _gammAfittog_equiyale^ 

weighting function becomes the one just as shown in equation (7). 



so 



A ft 

2 e, 2 = 2 { (x-f-b, y-t-c, z-M, ) M, ' } 2 -(7) 



55 



A part 4 .s a disturbance removal device for removing data errors in the weighted slice data from the 
weighting device 3 caused by disturbances such as the ones in which a secondary peak is generated at the 
tra.l.ng edge of a curve or the trailing edge of the curve is blunted as shown in Fig. 7 caused by the 
influence of blood ^-circulation etc. In Fig. 7. (,. 2. 3. - m -0 - n , are measurement data, and the data 
m + are the data of which trailing edge is blunted and data P are the data showing a secondary peak 
The curve connecting measurement data shows the curve of a gamma function to be found and the data m 



+ I, the data P, etc. are regarded to be the data which do not show original correct values because of 
disturbances. In the disturbance removal device 4, the percentage value RP for a peak value PV which 
appears in the measurement data is set beforehand. The measurement data m of the slice which fall 
smaller than (RP/100) x PV (peak value) for the first time will be found andnhe measurement data m are 
s made to be the final data to be used for gamma fitting. The disturbance removal device 4 transforms 
equation (7) to equation (8). 



*E= i { (x-t-bi y+ci z-Mi ) Mt ' } 2 -(8) 
to hi 

where m = final data to be used for gamma fitting 
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square value sum of residual to be evaluated 
A part 5 is a gamma function decision device for deciding a gamma function by finding the value of a 
variable which makes the square value sum E of the residuals in equation (8) obtained in the disturbance 
removal device 4 minimum. It performs operation as shown below. From equation (8). normal equations are 
obtained in which partial differential coefficients are put equal to zero as shown in the following. 



&E=0 
feE = 0 



These simultaneous equations are expressed as shown in equation (9). 

AX= B -(9) 

30 

Besides x. y and z. another parameter tO (refer to equation (3)) has to be decided to solve equation (9). For 
this purpose, in the device according to the present embodiment, the values of tO of several kinds to several 
tens of kinds, which are selected as explained later, are substituted for bl and d in equation (8). and 

35 equation (9) is solved independently for respective values of tO; in the result, a gamma function is decided 
from K, a, ^ and tO which make the square value sum E of the residuals minimum among all solutions for 
gamma fitting. How to select the several kinds to several tens of kinds of to values is shown in Fig. 8. A 
percentage value AP of measurement data for a peak value PV, and the number of times of repetition of 
fitting N are set beforehand. The measurement data (circle marks in the graph) of each slice are compared 

40 with the value (AP/100) x PV to find out a first slice of which measurement data become higher than 
(AP/100) x PV. and let T be the time of it. Let tO be a value selected in the following equation* and the 
gamma fitting of N times is performed. 
to = AT*j 
where AT = T.N. 

45 j = 0. 1. 2, ...N-1 

the result of gamma fitting of N times, a gamma function in which E becomes minimum is-decided-to- 

be the one to be found. 

The tO decided as described above is positioned as shown in Fig. 8. Let 1 be a slice right after tO, then 
the data of slices before I - 1 are not used. They are not used for gamma fitting, however, they are used for 
50 the evaluation of the square value sum E of the residuals as shown in equation (10) to be able to perform 
correct evaluation for any tO. 



Jt-r 

E= 2 M, 2 +■ £ { (x + b.y+c, z-M, ) M, ' } 2 

- (10) 



5 



EP 0 422 220 A1 



In equation (10). Ihe first term on the right side shows the slice data Which are not used for gamma 
fitting, and the second term shows the slice data which are used for gamma fitting and the slice data from i 
to m are reflected in equation (9). The gamma function decision device 5^ as described above solves a 
normal equation in equation (9) for tO's of N times, and a gamma function to be found is decided by K a b 
and tO when the square value sum E of the residuals become minimum. ' 

A gamma fitting device as described above has advantages as shown in the following 

(1) Calculation quantity is decreased much and calculation time is shortened much in performing least 
square approximation by solving normal equations obtained through logarithmic transformation of 
measurement quantities. 

(2) The gamma fitting equivalent to the gamma fitting in a real space can be performed by performina 
we.ght.ng .n a logarithmic space, which makes it possible to improve the accuracy of gamma fitting 

distunSTOe" 198 P0S$ib ' e t0 Perfarm aCCUrate 9arnma ' fittin9 by ori 9 inal data b V realizing the removal of 

(4) An optimum starting time of gamma fitting can be automatically decided by the decision method of to 
and the evaluation by equation (10). 

A gamma fitting device is not limited to the above-mentioned embodiment. In the embodiment 
explanation ,s given about the data, time and the increment of a CT number, however an arbitrary time 
sepal data can be approximated by a gamma function. As an example of a weighting function w OQ = X is 
adopted but another weighting function can be used. ' 

An image processor using the gamma fitting device as mentioned above is constituted as shown in Fig 
2 in Fig. 2, 1 1 is the gamma fitting device which performs gamma fitting operation for each picture element 
of the prctore element data of each tomograph image obtained in each slicing, and decides a gamma 
function wh.ch approximates most precisely to the time serial change of the picture data. A part 12 is a data 
control device which controls the gamma fitting device 11 to perform the operation for finding a gamma 
function fttted for respective picture element data based on the input data of scanned images; in the result 
all the sets of parameters K. a . » . and to which express gamma function curves of respective picture 
element data obtained from the gamma fitting device 11 are stored in the data control device by picture 
elements; and after that gamma function curves are reproduced by the above-mentioned parameters K a 
f' a " d ** various kinds of computer image data such as a peak value and an area below a qamma 

I" d6S f ^ ' ater - * rS ° Perated th9 ° peration pracess of «"Puter images as 
descnbed m the above, .f an excessive numerical value is generated and a picture element which may 
overflow.* found, the operation will be suspended and a mark of operation.suspension is stored in olace of 
the p.cture element data. A part 13 is an interpolation device, and when it finds a mark of operation 

SUJET? ? m M he inPUt d3ta fr ° m d3ta COntr0 ' d6ViCe 12 ' re " Search of ima 9 e data is P erf °™ed over the 
range of N x N matnx centering the picture element, and when there are found more than M effective data 

in the matnx of N x N. an average value of these effective data will be substituted for the mark of operation 

suspension to fill data blank, and when picture element data are completed, a computer image constituted 

by the data control device 12 is supplemented to be completed. By this arrangement, it is made possible to 

supplement smoothly the value of an incalculable picture element caused by overflow, etc. by interpolation 

from the crcumference and to obtain a computer image of good picture quality and good to look at A 

computer .mage output from the interpolation device is displayed in a display device (not shown in ihe 

drawing . The matrix size N and the number of data M shall be decided properly depending on the olace of 

observation or the degree of a morbid state. The data control device 12 calculates various kinds of 

computer .mages, for example, as shown in the following which are specified beforehand based on the 

gamma fitting performed by the gamma fitting device 1 1 . 

— (t)-A-peak-value-of-a-garnma-functioTi curve— — - - - - 

(2) A period of time from the origin of the coordinate to the peak value 

(3) An area below a curve 

(4) Inflection point width 

(5) Time of appearance 

(6) Mean passage time 

(7) Corrected mean passage time 

(8) Relative flow rate 

The operation is suspended for a picture element of which overflow is brought about in the orocess of 
above-mentioned operation of computer images, and a mark expressing the suspension is giveHo the 
P-ctore element. When the interpo.ation device 13 finds the mark, i, tries to obtain an average value of more 
man m effective data out of N x N data in the vicinity of the mark for interpolation. anS c omplete 7Z 
operation of a computer image of which the operation is not completed including the interpolated data to 



display the image in a display device. 

The explanation on a most preferred form for the execution of the present invention is given as 
described in the above, however it will be easy to make various transformation without departing from the 
range of the following claims for the persons having ordinary knowledge in the technical field to which the 
s present invention belongs. 



Claims 



to 1. An image processor provided with a gamma fitting means for setting a gamma function curve 
approximating to a curve showing the movement of a contrast medium in a tomogram image obtained 
by a tomograph and for constituting various characteristics of the set gamma function curve as a 
computer image comprising: a data clip means for performing coordinate transformation making the 
leading data of measurement of an image changing with the lapse of time accompanying the injection 
is of a contrast medium be the origin of the coordinate, and for shifting the data being made to negative 
data by the coordinate transformation to positive data; a logarithmic transformation means for storing a 
3 element first degree equation, a gamma function undergone logarithmic transformation, and for 
substituting the data from said data clip means for the coefficients: a weighting means for weighting the 
output data of the logarithmic transformation means; a disturbance removal means for removing 
w disturbance data to correct the disorder of a curve caused by disturbances; and a function decision 
means for setting an approximate gamma function curve by obtaining unknown parameters of a gamma 
function by operation based on the data after removal of disturbance from said disturbance removal 
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means. 



An image processor for setting a gamma function curve ap proximating to a curve showing the 
movement of a contrast medium in a tomogram image obtained by a tomograph, and for constituting 
various characteristics of the set gamma function curve to a computer image comprising: a data clip 
means for performing coordinate transformation making the leading data of measurement of an image 
changing with the lapse of time accompanying the injection of a contrast medium be the origin of the 
coordinate, and for shifting the data made to negative data by coordinate transformation to positive 
data: a logarithmic transformation means for storing a 3 element first degree equation, a gamma 
function undergone logarithmic transformation, and for substituting the data from said data clip means 
for the coefficients: a weighting means for weighting the output data of the logarithmic transformation 
means: a disturbance removal means for removing disturbance data to correct the disorder of a curve 
caused by disturbances; a gamma fitting means comprising a gamma function decision means for 
setting an approximate gamma function curve by obtaining unknown parameters of a gamma function 
by operation based on the data after removal of disturbances from said disturbance removal means: a 
data control means for forming a computer image based on a gamma function curve obtained by the 
calculation of gamma function decision parameters output from said fitting means by controlling the 
operation of said gamma fitting means;and an interpolation means for completing a computer image by 
interpolating the picture elements in an incalculable part in the operation process by said data control 
means. 
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